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ABSTRACT 

Ice crystal clouds in the upper troposphere can generate polarisation signals at the 
|j.K level. This signal can seriously affect very sensitive ground based searches for E- 
and B-mode of Cosmic Microwave Background polarisation. In this paper we estimate 
this effect within the CLOVER experiment observing bands (97, 150 and 220 GHz) for 
the selected observing site (Llano de Chajnantor, Atacama desert, Chile). The results 
show that the polarisation signal from the clouds can be of the order of or even bigger 
than the CMB expected polarisation. Climatological data suggest that this signal is 
fairly constant over the whole year in Antarctica. On the other hand the stronger 
seasonal variability in Atacama allows for a 50% of clean observations during the dry 
season. 

Key words: Cosmic microwave background, cosmology: observations, techniques: 
polarimetric, atmospheric effects 



1 INTRODUCTION 

CLOVER (C^ObserVER) is a collaboration between the 
Cardiff Astronomy Instrumentation Group, Oxford Astro- 
physics, Manchester Astrophysics and the Cavendish As- 
trophysics Group in Cambridge, on an experiment to mea- 
sure the Cosmic Microwave Background (CMB) polarisation 
( Maffei et aTp004 l. 

Polarisation of the CMB is caused by Thomson scat- 
tering of CMB photons at the last scattering surface (|Hu| 



& White 19971. The signal can be decomposed into a curl 



and a curl-free component, known as B- and E-mode. The B- 
mode signal, which is at best one order of magnitude weaker 
than the E-mode, is generated by primordial tensor pertur- 
bations and therefore its detection would provide valuable 
information about the history of the early universe. 

The main ambitious objective of CLOVER is the mea- 
surement of the B-mode; in order to achieve this result the 
experiment will deploy large format imaging arrays, oper- 
ating at 97, 150 and 220 GHz with 30% bandwidth and a 
beamwidth of approximately 8 arcmin; the instrument is de- 
signed with an unprecedented level of systematic control and 
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will be deployed in the Atacama desert (Chile) at an alti- 
tude of 5080 m. Alternative sites are the Antarctic stations 
of Dome C and South Pole. 

In spite of a site choice with favorable atmospheric con- 
ditions we expect the signal of the atmospheric fluctuations 
to be well above the intrinsic instrumental noise. 

Since the CLOVER receiver modulates signal polarisa- 
tion, the main concern about atmospheric effects is about a 
potentially polarised signal from the atmosphere. Water va- 
por is the major absorbing component at mm wavelengths 
and its spatial distribution is highly variable with time. 

These variations could also introduce some polarisation 
noise; in situ measurements of the turbulence suggest that 
this polarised contribution to system noise is expected to be 
gaussian and negligible during most of the observing time, 
even in Atacama which should be the worst of the three 
sites, both because of the stronger day-night thermal cycle 
and the height of the mean boundary layer that can be from 



200 to 2000 meters (Giovannelli et al. 20011, with respect 



to 230 m in South Pole and 30 m in Dome C (Agabi et al 



2006 1 during winter. 



In addition to the variable contribution by water vapor, 
the strong oxygen features at 120 GHz and around 60 GHz 
dominate the brightness temperature of the atmosphere in 
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Ce OVER's spectral region. The presence of the Earth's mag- 
netic field causes a Zeeman splitting of the energy levels of 
the oxygen, thus resulting in a polarised emission depend- 
ing on the relative alignment between the line of sight and 
the magnetic field. This effect is well known for atmospheric 

and references 



Table 1. Candidate CLOVER observation sites. 



measurements (see von Engeln et al 



cited therein). Keating et al. (19981 discussed its impact 



(1998), 



on CMB measurements. Hanany & Rosenkranz (20031 esti- 



mated that the circularly-polarised component is not neg- 
ligible if the intrinsic leakage between linear and circular 
polarisation in the instrument will be of the order of one 
per cent. 

However, the polarised intensity due to oxygen is not 
expected to vary with time but is be fixed for a particular 
azimuth and elevation direction. Hence any scanning strat- 
egy will modulate any residual atmospheric signal in a very 
predictable way. Also, the estimates are for the DC level of 
the signal. Oxygen is well mixed in the atmosphere at alti- 
tudes up to approximately 80 km, hence fluctuations in the 
oxygen signal on the angular scales to which CLOVER is 
sensitive will be very small. It should therefore be possible 
to separate this signal from the CMB polarisation well down 
below the sensitivity required. 

Upper tropospheric ice clouds (like cirrus clouds) repre- 
sent another source of polarised radiation. These clouds are 
at high altitudes and contribute to the energy budget of the 
atmosphere (greenhouse effect) since they absorb thermal 
IR radiation from the ground and, as they are cold, emit lit- 
tle infrared radiation. This warms up the Earth-atmosphere 
system. On the other hand, ice clouds reflect incoming solar 
short wave radiation and hence cool the Earth-atmosphere 
system. 

At mm and sub-mm wavelengths the interaction be- 
tween ice clouds and radiation is mainly due to scattering. 
Absorption is negligible, and so is the thermal emission. The 
scattering by ice clouds will introduce a polarisation sig- 



nal. Teichmann et al. ( 20061 have shown that this polarisa- 



tion signal arises even assuming spherical ice particles, due 
to the asymmetry of the radiation field in the atmosphere. 



However, real cloud ice particles are not spherical (Wallace, 



& Hobbs 19771, and this increases the polarisation signal. 
Moreover, there is a growing evidence of horizontal align- 
ment of cloud ice particles due to a combination of aerody- 
namic and gravitational forces ( Prigent et al.||2005 l, which 
further increases the polarisation signal. The actual magni- 
tude of the cloud polarisation signal will depend strongly on 
the particle size and shape, and on the line of sight direction. 

Experimentally, ice crystal depolarisation is a well- 
known problem for high frequencies satellite telecommuni- 
cations (20-^50 GHz) based on signal polarisation diversity 
encoding (Martcllucci et al. 20021. Measurements carried 
out with experimental telecommunication payloads (ITAL- 
SAT and OLYMPUS) showed that even at relatively low 
frequencies (with respect to CLOVER bands) the depolaris- 
ing effect of ice crystals is not negligible ( Trione|2 003). 

For CMB measurements, the impact of the cloud scat- 
tering is twofold. Firstly, the CMB signal is depolarised, sim- 
ilarly to a telecommunication signal. Secondly, the cloud 
also scatters back upwelling thermal radiation from the 
earth surface into the line of sight of the instrument. For 
telecommunication links this second effect is negligible, due 
to the large intensity of the telecommunication signal. But 



Site name 


Latitude (°) 


Longitude (°) 


Height (m) 


Atacama 


23 S 


67 W 


5080 


Dome C 


75 S 


123 E 


3280 


South Pole 


90 S 




2900 



for CMB measurements the radiation scattered back by the 
atmosphere will often be more intense than the CMB signal. 
The backscattered signal will be partially polarised, and its 
polarisation characteristics will depend on many factors, as 
will be explained in the following sections. 



2 ICE IN THE UPPER TROPOSPHERE 
FOR THE THREE SITES 

Although ice clouds play an important role in the atmo- 
sphere energy budget, up to now they are poorly measured 
and modeled. 

Satellites provide global measurements of integrated ice 
mass (Ice Water Path - IWP) with frequent revisit time on 
a long term basis. Sensors detect both reflected sunlight (in 
the UV and visible, Buri ez et al.| ( [T997| ) and thermal emis- 
sion. While the last method is limited to semitransparent 
clouds the first one only works if the earth surface albedo 
is not too high (which, unfortunately, is exactly the case for 
Atacama desert and Antarctica). 

The MODIS instrument carried by the Terra and Aqua 
EOS-NASA satellites is equipped with a cirrus clouds detec- 
tion band at 1.38 um; the method, first suggested by |Gao et] 
al. ( 1993 1 for airborne measurements, suffers the same draw- 



back since it assumes that upwelling radiation reflected by 
the earth surface is strongly absorbed by water vapor in the 
lower troposphere and therefore the method is not effective 
when water vapor column density is very low (i.e. < 4 kg/m 2 
corresponding to 4 mm precipitable water vapor, PWV). 

The next generation satellite or air-borne instrument 
will characterize ice clouds by measuring from above the 
radiation brightness temperature depression with imaging 
radiometers in the mm and sub-mm range, as proposed 
by |Buehler et~aL] (2005a} and |Evans et at] (|2002|| and 



will include also polarisation measurements (Hayton et al 



20031. In situ and aircraft-borne experiments provide the 
most accurate information on ice clouds. Dual polarisation 
radars (30-90 GHz), polarisation diversity lidars and air- 
planes equipped with cameras are used to characterize ice 
density, crystal shapes, orientation and size distribution. 

As mentioned above, there is evidence that ice needles 
and plates (especially those with large size and aspect ra- 
tio) have a preferred orientation. As crystals drift down- 
wards, they become oriented in a maximum drag condi- 
tion: aerodynamic forces tend to cause their long axis (or 
axes) to fall horizontally (i.e. the shortest axis is perpendic- 



ular to the ground; see Evans et al. ( 1998 1 and references 



therein for both models and some experimental results). 
Such a good orientation is testified also by the relatively fre- 
quent presence of optical effects, such as sun haloes, which 
happens only if the crystal are aligned within few degrees. 
More recently, |Noel fc Sassen| (2005 ) derived this result from 
polarised lidar backscatter measurements. |Noel fc Chepfer] 
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Figure 1. IWP percentiles (values express column density, units are g/m 2 ); each graph represents a three month period. IWP values 
for the 25% percentile in Atacama, during JJA and SON, are $J 0.0001 g/m 2 . 



(20041 analysed polarised visible light measurements from 



the POLDER satellite, and found that 50% of high clouds 
show a glint signature implying at least a fraction of the ice 
platelets to be horizontally aligned (within a very narrow an- 
gle). [Dwis^etjir (20051 showed that cloud ice crystals gen- 
erate a polarisation signal in the limb measurements at 122 
GHz (carried out by the Microwave Limb Sounder (MLS) 
on the Aura satellite). They also concluded that the effec- 
tive particle shape can be approximated by a horizontally 
oriented oblate spheroid with an aspect ratio of 1.3. This 
rather moderate value of the effective asphericity is due to 
the fact that, while there are individual particles with large 
aspect ratios, there is also an averaging effect over the dif- 
ferent sizes, shapes and orientations of the individual ice 
particles. 

The purpose of this paper is to estimate the influence 
of cloud ice particles on the CLOVER measurements at the 
selected and at the two backup sites (Table 1). For this pur- 
pose assumptions are made on the range of cloud ice amount 
to be expected for the different sites, as well as on the ice 
particles size, shape, and orientation. 

We used the general circulation model of the European 
Centre for Medium Range Weather Forecasting (ECMWF, 



see Uppala et al. (2005 1 ) to estimate the statistics of cloud 



ice mass. Figure 1 shows the statistics in the form of per- 
centiles. It shows for example that Atacama, though gener- 
ally a very dry place, can have ice clouds exceeding an Ice 
Water Path (i.e. ice column density, IWP hereafter) of 100 
g/m 2 for about 10% of the time from December to May. 
These statistics were derived from a three years long (2000 
to 2002) data set obtained with a gridded version of the 
ECMWF model data with a grid resolution of 1.5° by 1.5°. 
A longer data set is still under analysis, however we believe 
that the period considered is quite significant since it was 
not affected by climatic extremes (such as El Nio events in 
South America). 

As expected, Atacama shows a much stronger seasonal 
variability than Antarctica, in the driest season IWP is of- 
ten (25% of time) ^ 0.0001 g/m 2 . Since this is the selected 
site for CfOVER experiment, we evaluated the polarised 
signal from ice crystals using a standard atmospheric pro- 
file for these latitudes, re-scaled in accordance with locally 
measured water vapor values (1 mm PWV, i.e. 1 kg/m 2 ). 

The impact on CMB polarisation measurements will be 
discussed for a wide range of IWP values from 0.0001 g/m 2 
to 100 g/m 2 . 
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Figure 2. Simulated cloudy radiances, radiance difference and 
polarisation difference between a cloudy and a clear sky as a func- 
tion of instrument zenith angle. Different curves are for different 
observation frequencies. 



3 SCATTERING MODEL 

There are no direct data available for particle size, shape 
and orientation. A realistic size distribution based on litera- 



Ivanova et al.[2001 


and 


1 g/m^ (Donovan 


2003 



All 



columns, for IWP < 1 g, 
pact polycrystals for IWP ^ 
particles were horizontally aligned with random azimuthal 
orientation. The equal-volume-ellipsoid aspect ratio of the 
particles is assumed to be 1.3. 

The radiative transfer (RT) model used for this study 
is able to handle also even more realistic cases, but for this 
first assessment it was decided to keep the assumptions as 
simple as possible for clarity. 

The RT model used was the Atmospheric Radiative 
Transfer Simulator (ARTS). The basics of the model are 
described by |Buehler et al.| ( |2005a[ ) . Here we used version 
ARTS-1-1-1095, which can simulate the scattering of radia- 
tion by cloud particles ( |Eriksson et al.||2005 l. ARTS offers 
two different scattering algorithms: a Monte Carlo algorithm 
and an iterative discrete ordinate algorithm (DOIT). In this 
work we used the DOIT algorithm, which is described in 
detail by |Emde et al.] p004| . 

The clear-sky part of ARTS has been compared against 
a range of other microwave radiative transfer models 
1 Melsheimer||2005 1 and against co-located AMSU data and 
radiosonde profiles ( |Buehler et al.|2004| ). The scattering part 
of ARTS has been compared against several other scatter- 



ing models (Emde 2005 Hoepfner & Emde 20051, against 



co-located AMSU data and mesoscale weather prediction 
model fields ( Sreerekha et al.|2005 1. 



Figure 3. Simulated radiance, radiance difference and polarisa- 
tion difference between a cloudy and a clear sky at a zenith angle 
of 5° as a function of IWP. 



ARTS can handle all four Stokes components. However, 
the azimuthally symmetric geometry in this case implies 
that only I and Q Stokes components are non-zero. The 
component I represents the total intensity (sum of horizon- 
tally and vertically polarised intensities), the component Q 
represents the linear polarisation difference. 



4 RESULTS 

Figure 2 shows the ARTS simulation results. The sensor was 
assumed to be at an altitude of 5080 m. The cloud was as- 
sumed to be located at an altitude of 9000-11000 m. Each of 
the plots includes three curves corresponding to the central 
frequencies of the CLOVER bands. The assumed IWP value 
is IWP = 1 g/m 2 . The different rows show the simulated 
radiance (top row), the radiance difference (i.e. difference 
between the radiance from the cloud and the clear sky, mid- 
dle row) and the polarisation difference (Stokes component 
Q, bottom row). All signals are shown as a function of the 
instrument angle relative to the zenith direction. The polar- 
isation signal is zero at the zenith, and increases with zenith 
angle. 

Both the radiance difference and the polarisation signal 
increase with the IWP. This is shown more clearly in Fig- 
ure 3, where the simulated radiance, the radiance difference 
and the polarisation difference are plotted as a function of 
IWP for a zenith angle of 5°. The figure covers values of 
IWP ranging from 0.01 to 100 g/m 2 . Note that the radia- 
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Table 2. Simulated polarisation signal (Q) per unit IWP at three different zenith angles. For every angle the total signal and the 
polarisation induced on 2.7 K CMB are quoted. 
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Figure 4. Angular power spectrum of CMB polarisation (E-mode and B-mode) as calculated by CMBFAST code. Straight lines represent 
an upper limit on the amount of polarization induced by ice crystal clouds on the 2.7 K CMB for more than 50% observing time during 
dry season in Atacama (IWP = 0.01 g/m 2 ) at the three CLOVER frequencies and an upper limit for 25% observing time during dry 
season (IWP = 0.0001 g/m 2 ) at 97 GHz. A flat power spectrum is assumed; although it appears that the ice crystal signal might be 
dominating ground-based observations of CMB polarization of E and B-modes, it must be stressed that, while the CMB signal is fixed 
in the sky, the ice signal is most probably variable with time. Therefore it is always possible to disentangle and therefore greatly reduce 
the ice signal from the sky signal by a properly designed observing strategy 



tive transfer is in the linear regime at these frequencies for 
reasonable IWP values. (Note the logarithmic x-scale of the 
plots) . 



5 EFFECT ON OBSERVED CMB 
POLARISATION 

The results of the calculations using ARTS show the ex- 
pected polarisation signal induced by the ice clouds. The 
signal represents the difference between a cloudy and clear 
sky. Cloud parameters are derived from climatological data. 

This signal is divided into two components: 



© 2002 RAS, MNRAS 000,[T}{7] 



6 Luca Pietranera et al. 



(i) upwelling radiation, emitted by ground and lower at- 
mosphere, that is back-scattered and polarised by the ice 
crystal clouds; 

(ii) CMB radiation that is forward-scattered by the ice 
with a variation in its polarisation (ARTS model assumes 
that CMB radiation is unpolarised). 

In order to disentangle the two components, we run the 
ARTS model alternately imposing T 'ground and Tcmb equal 
to °K. Table 2 shows the polarisation signal per unit IWP, 
which is applicable in the linear regime (at least up to 100 
g/m 2 )- 

In order to compare the ice crystal induced polariza- 
tion with the expected cosmological signal, we generated the 
angular power spectrum of the CMB polarisation (both El- 
mode and B-mode) in an ACDM (Adiabatic Cold Dark Mat- 
ter) cosmology scenario with tensor to scalar ratio r = 0.01 



(using CMBFAST code, see Seljak fc Zaldarriaga (19961) 



Figure 4 reports the expected CMB polarisation signal 
power spectrum in terms of the spherical harmonic coeffi- 
cient i. The range is limited to spatial scales which are rel- 
evant for CLOVER experiment: between 10° and 10 arcmin 
(£ ranging from 20 to 1000). 

The polarisation signal on CMB radiation due to ice 
crystals is plotted for each CLOVER frequency, assuming an 
observing zenith angle of 5° and IWP values corresponding 
to an upper limit for respectively 50% and 25% of observing 
time during dry season. 

In the figure we arbitrarily assumed a flat power spec- 
trum for the observed signal from ice crystal clouds; it has 
to be underlined that this model is really inadequate and 
a more realistic representation should take into account the 
cirrus cloud morphology and spatial distribution which un- 
fortunately are very poorly known. Qualitatively, a decrease 
at the high £ is expected, since the clouds are supposed to 
be quite homogeneous at these angular scales, but, in any 
case, the signal could still be orders of magnitude stronger 
that the expected B-mode component during a significant 
fraction of the observing time. 



6 CONCLUSIONS 

The polarisation signal on CMB radiation due to ice crystal 
clouds cannot be neglected by a ground based experiment 
looking for mapping E- and B-mode patterns even for low 
ice water column density values (0.001 g/m 2 ). The effect is 
particularly strong for the high frequencies commonly used 
for CMB measurements (150 GHz and above). 
Possible mitigation strategies include: 

(i) constant elevation fast instrument scanning: at the 
moment CLOVER is designed to internally modulate the 
CMB polarisation signal and therefore an absolute measure- 
ment of this parameter (i.e. for each pixel) would be possi- 
ble. A differential approach, based on rapidly measurement 
of signal differences between contiguous pixels would can- 
cel the constant bias due to ice crystals clouds. However, 
the measurement would still be affected by the cloud spatial 
distribution and inhomogeneities; 

(ii) Ice Water Path could be measured by processing 
Earth Observing (EO) satellite data. In principle geo- 
stationary satellites near real-time monitoring capabilities 



(one measurement every 15-30 minutes). However, as dis- 
cussed above (para. 2), very low tropospheric water vapor 
content (both in Atacama and in Antarctica) seriously un- 
dermine this measurement. Future EO mission dedicated to 
ice cloud measurements will be helpful for solving this prob- 
lem; 

(iii) IWP in situ measurements can be carried out in or- 
der to assess the quality of measured polarisation, and to 
characterize, during clean nights, a set of reference pixels 
to be used for ice detection. Among the possible options, 
the use of a polarised lidar, though difficult to deploy and 
operate, would provide us with a full characterization of ice 
crystals (shape, size distribution and orientation) allowing 
an accurate modelling of radiative effects in the microwaves. 

In this difficult observational context, a preliminary 
comparison between the sites can be attempted. The pres- 
ence of the sun at high elevation angles penalizes Atacama 
with respect to Antarctica sites; however, due to the high 
sensitivities required for CMB polarisation measurements, 
it is in any case difficult to carry out measurements with 
the sun above the horizon (strong signal in the instrument 
sidelobes). On the other hand climatological data on ice 
cloud occurrence and density suggest that Atacama observ- 
ing conditions during dry seasons are significantly better 
than Antarctica (with frequent occurrences of IWP ^ 0.0001 
g/m 2 ) 
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